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Abstract 

^^ . We present a program to calculate the total cross section for top-quark pair production 

H I in hadronic collisions. The program takes into account recent theoretical developments 

such as approximate next-to-next-to-leading order perturbative QCD corrections and it 
allows for studies of the theoretical uncertainty by separate variations of the factoriza- 
tion and renormalization scales. In addition it offers the possibility to obtain the cross 
section as a function of the running top-quark mass. The program can also be applied 
to a hypothetical fourth quark family provided the QCD couplings are standard. 
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1 Introduction 

The top-quark is the heaviest elementary particle in nature discovered so far. As a consequence 
of its large mass close to the scale of the electroweak symmetry breaking it has remarkable 
properties making it a distinct research object. For example the short lifetime does not allow 
the formation of hadronic bound-states. Rather, the top-quark decays before it hadronizes, a fact 
often referred to colloquially as the top-quark behaving like a quasi-free quark. Non-perturbative 
effects are thus essentially cut-off by the short lifetime and, as an important consequence, the 
polarization of top-quarks can be studied through the parity violating decay into a W^ -boson and 
a bottom quark. The top-quark also provides an interesting environment for precision tests of the 
Standard Model (SM) and possible extensions, e.g., by constraining the allowed range for the 
Higgs mass. 

A mandatory ingredient for top-quark physics at hadron colliders are precise theoretical pre- 
dictions to compare with. Current Tevatron measurements and the perspectives at LHC, i.e., a 
measurement of the top-quark pair cross section with an uncertainty of the order of 5% only set 
the target for theoretical predictions of the production process. Clearly such an accuracy needs 
to include quantum corrections. 

Within Quantum Chromodynamics (QCD) radiative corrections were calculated some time ago 
to next-to-leading order (NLO) first considering unpolarized top-quark production Id 13 and 
later including spin information Q . In the former case, also completely analytical results have 
recently been provided DUl. Beyond the NLO accuracy in QCD various sources of possible im- 
provements have been identified. Large logarithmic corrections due to soft gluon emission were 
investigated and resummed at the next-to-leading-logarithmic (NLL) accuracy (SMl and recently 
improved to include also the next-to-next-to-leading-logarithmic (NNLL) corrections JTHll. Al- 
ternatively, resummation has also provided the means to construct parts of the full next-to-next- 
to-leading (NNLO) fixed order results [|71ll0l[IT||, which can be supplemented by including all 
Coulomb type corrections IfTTI and also the full scale dependence at NNLO accuracy ir71[T2]|. 
With a target precision for the total cross section at the few per cent level also bound state effects 
from the resummation of Coulomb type corrections lfT3HT4| as well as electro-weak radiative 
corrections at NLO ifTSUTTl need to be considered. 

The compilation of all these results is in principle straight forward given the extensive literature 
on the subject. However, no publicly available program exists so far which contains the latest 
theoretical developments. The "modus operandi" of the past was that predictions were updated 
by theorists from time to time taking into account new theoretical improvements and/or new sets 
of parton distribution functions (PDFs). The aim of the present paper is to provide a program 
for the computation of the top-quark pair cross section including state of the art theory. As such 
it can serve as a reference for future cross section calculations. The program Hathor includes 
perturbative QCD corrections at higher orders in the different approximations along with options 
allowing also a detailed study of the theoretical uncertainties. Moreover, it provides the possi- 
bility to compute the total cross section not only in the commonly adapted pole mass scheme 
but also in terms of the MS mass a choice recently employed for top-quark pair production in 
hadronic collisions for the first time lfT2l[T8ll . Finally, the aim of this publication is not only to 
provide a tool for cross section calculations but also to facilitate experimental analyses. To that 



end the package contains in addition to the stand alone program also a small library that can be 
easily integrated into existing code. 

The outline of the article is a follows. In the next Section we briefly discuss the theoretical 
foundations as well as the procedure to convert to the MS mass. Sections [3] and |4] contain instal- 
lation details and the program description while usage and examples are given in Section [5l We 
end with conclusions in Section [6l All formulae as implemented in Hathor are collected in the 
Appendices |A] and |Bj 

2 Methods 

The hadronic cross section for top-quark pair production is obtained from the convolution of the 
factorized partonic cross section d;/ with the parton luminosities Lty. 

s 
<^h^h2^ttx{S,mt) =Y^ ds Lij{s,S,Hf) Oij{s,mt,as{nr),l^f) , (1) 

4mf 

Lij{s, S,fJf) = -J jfi/h^ ( i'^^ ) ^J/''^ (?^V ■ ^^^ 

s 

Here S denotes the hadronic center-of-mass energy squared, /Jr, (jj/) denotes the renormalization 
(factorization) scale and the functions /,y;,j ^ i^il-Jf) are the PDFs describing the probability to find 
a parton of type i with a momentum fraction between x and x + dx in the hadron h^. The QCD 
coupling constant UsiiUr) is evaluated at the scale fXr- In the following we use a^. in the scheme 
with Hf light flavors. For top-quark production, the running is thus determined by the five light 
flavors u,d,c,s,b which we treat as massless. The top-quark mass rut appearing in Eq. ^ is the 
mass renormalized in the on-shell (pole-mass) scheme. 

In perturbative QCD the partonic cross section Oij{s, mt, (ls-,Hf) is expanded in the QCD coupling 
constant up to NNLO: 

^ij = a^ ^ij (s^rnt) + al 6}j' {s,mt,fjr,fjf) +a'l6}f{s,mt,iJ,;l^f) + O(a^) , (3) 

with as = cx.v/7r. 

In leading-order (LO) only the parton channels qq and gg contribute and the respective Bom 

cross sections are given by: 

*£' = |^>(3-P^), (4) 

*» = |^^{(33-I8|5^ + P*)l„(i±|)-59|5 + 31P'}, (5) 

with P = y/l — p and p = 4mj/s. Starting from NLO also the gq and gq channels contribute. In 
Ref. fH (and in many subsequent publications) an alternative decomposition was used in terms 



of so-called scaling functions fif. 

2 

(6) 
Since the scaling functions are dimensionless they depend only on p and the ratios fjf/nit and 
/Jr/lJf- The full renormalization and factorization scheme dependence can be constructed using 
the renormalization group equation, the standard evolution equations of the PDFs and informa- 
tion about lower orders, i.e., up to NNLO knowledge of /, ■■ (p) and /j- ■ (p, 1, 1) is sufficient. 
The general structure can be written in the following form 

4^^(p,A//A%,A/.M) = fll'^+LMf!j'^ + 20Lj,fj^, (7) 

+2lLRfj^^ +30^Llfjf\ (8) 

with ij = {qq, gg} and we abbreviate Lm = \n{/j^Jmj) and Lr = \n{fi-^/fA). The scale dependence 

in the gq (gq) channel can be easily derived from the above realizing that fgq' = so some terms 
in Eqs. (|7]) and dS]) are simplify absent. In the conventions used here the coefficients of the 
beta-function 0, 1 are given by 

The Bom contributions have been presented in Eqs. ©, ([5]) and at present also the complete 
NLO corrections are known, i.e., the functions /^ ■ and //■ in Eq. ([7]). A complete NNLO 

calculation for the total cross section is not yet available, since f^j '{p,/Uf/mt,iLir/inr) in Eq. ^ 

(^0) (21) (22) 

is missing the contribution />" ^ while /^ and /) ■ ' have been obtained from renormalization 

group arguments as mentioned above. However, exact expressions for f-^- ' in the limit p — >> 1 
based on soft-gluon resummation have been derived and provide the foundation for approximate 
NNLO results o{<5hih2^ttx- 
The central physics questions to be addressed can be phrased as follows: 

• How large is the total cross section <5h^h2^ttx at a given order in perturbation theory ? 

• Given a computation of the total cross section according to Eq. ([T]) what is the associated 
theoretical uncertainty ? 

In order to address these issues the package Hathor has different production models implemented 
which are accessible to the user as options. In the following we briefly describe these options 
as far as the underlying physics is concerned. To be self-consistent and for easier reference all 
necessary theory input, e.g., the scaling functions has been collected in Appendix [Hi For details 
of how to access these options when running Hathor we refer to the next Sections H] and [51 



Option LO 

The option LO provides a rough estimate although with large theoretical uncertainties which 
will receive sizable corrections at higher orders. This option uses the Born cross sections of 
Eqs. dH), © (see also Eqs. (IB11-(IB3])). 



Option NLO 

The option NLO is the first instance where a meaningful theoretical uncertainty can be quoted in 
perturbation theory. This option employs the complete NLO QCD corrections [[Bill. All scaling 
functions /-^ ' are given as accurate fits lfT2ll based on the recently published analytic results (H, 
(see also Eqs. (iRih-dRel)). 



Option NNLO 

The option NNLO is required whenever predictions with an uncertainty of better than 0(10)% are 
needed. This option is based on the known threshold enhancement due to soft gluon emission, 
i.e., complete tower of Sudakov logarithms at NNLO accuracy, supplemented by all Coulomb 
type corrections IfTTI and also the full scale dependence |I71[l2l (see Eqs. (IB.12I) - (IB.17I) ). This 
ansatz provides a good approximation for the total cross section [TTllTOl. a fact which is supported 
by the observation that the QCD corrections to top-quark pair production in association with an 
additional jet are small lfT9l - BT]| . 
Thus, using the results of Refs. [|11[I12[| . we have for the functions /,■ : 

4? = T^^{^ln4|3 + (-1505.1589 + 37.925926n/)ln3p (10) 

+ (1046.4831-90.838135^/- 140.36771^) In^ (3 
+ (249.67547 + 55.776275n/ + (54.038454 - 4.3864908/7/) -) In p 

+3.6077441-2 + (-5.2728242 + 1 .8447758n/) ^ + C^J } , 

.(20) _ _P^657r 3 2n (... 

^^^ - (16^^^"^^'^^' ^^^^ 

A? = 7T^^{46081n4|3+ (-2321.5810 + 85. 333333n/)ln3p (12) 

+ (-315.57218- 119.35529n/ + 496.30011-)ln2p 
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+ (2346.8889 + 21.969529n/+(286.67132 + 6.8930570n/)^jlnp 



+68.547138p + (-3.7910584 -0.96631115n/)^ +4?}, 

where the unknown functions Q= and Cgg in Eqs. (flOl) and (fT2l) parametrize the contributions 

which are not enhanced in the threshold region, i.e., 0{^^). The g^-channel, that is fgq in 
Eq. (fTTI) . is additionally suppressed near threshold with corrections of order 0(|3^ In^(P)) . 
In summary, the option NNLO (which has been used e.g., for the phenomenological studies of 
Ref. (T2^) uses all presently available information at NNLO. In this way, it attempts to construct 
the relevant parts of the complete NNLO corrections. Necessarily, the small associated theoret- 
ical uncertainty [J71[l0l| due to scale variation (jjr and /jf) estimates effects beyond NNLO. An 
additional systematical uncertainty on the quality of the approximate NNLO result can be quan- 

(2) (2) 

tilled by varying the constants Q= and Cgg in a reasonable range comparable to the size of the 
other coefficients in Eqs. ^ and ^, i.e., C^f = ±0(100). The defauh value is C^f = 0. 



Option LOG_ONLY 

The option LOG_ONLY is also motivated by the idea of soft gluon enhancement near threshold. It 
emerged as a conservative definition of the theoretical uncertainty in a comparison of different 
approaches to incorporate dominant terms beyond NLO to a certain logarithmic accuracy. In 
particular threshold resummation at NLL accuracy, performed as in Refs. [I6ll22l| which typically 
proceeds in Mellin-space, see Eq. (fT3] ). has been tested against an expansion in powers of ln^(|3) 
in momentum space as advocated in Ref. [TTlfTOll. This comparison has yielded satisfactory agree- 
ment, because resummation beyond NLL, i.e., at NNLL accuracy has only a minor effect |[7|. 
The option LOG_ONLY as discussed below is based on work with CCMMMNU |l23l|. It is a gen- 
uine NLO approach with logarithmic improvement near threshold and scale variations in /Jr and 
/Jf (with a constraint on the ratio of /Jrf/Jf) estimate effects beyond NLO. Being a conservative 
approach the resulting theoretical uncertainty is necessarily larger than in the option NNLO. 
Let us briefly mention the essential technical points. In Refs. [|6ll22l the logarithmic enhancement 
is constructed from the ln(A'^) terms in Mellin space where the resummation is usually performed. 
The transformation from momentum (p— ) space to Mellin (A'^-) space is given by 

a{N)= f Jpp^-ia(p), (13) 

Jo 

where p = 4mj/s. The important feature of Eq. (fT3l) to realize is that beyond logarithmic ac- 
curacy the momentum space and the Mellin-space expressions do differ by terms which are not 
enhanced in |3 or, respectively power-suppressed in l/N. Any difference could be included in a 

choice of C^f in Eqs. ([H), ^. 

Option LOG_ONLY has to be used with the Option NNLO. It applies Eqs. ([JO]) and (O, but truncates 

the function /■■ beyond NLO to is leading term in (3 (cf. Eqs. (H)), ([5])) and, for consistency 

neglects the g^-channel beyond NLO (fgq and Eqs. (IB. 141 ) and (IB. 151 )). Also for the scale 
dependent part, the option LOG_ONLY keeps only terms which are logarithmically enhanced in the 



threshold region. In this case the functions /^ and /j ■ in Eqs. (|B.12I) - (IB.17I) ) are truncated 

(2) (2) (2) 

to logarithmic accuracy. Again, one could also vary the constants C g and Cgg in a range Q = 

(2) 

±(9(100) to test for additional systematical uncertainties. The default value is C- — 0. 



Option MS_MASS 

The option MS_MASS allows for the computation of the total cross sections as a functions of the 
running mass in the MS scheme. In a nut-shell this is based on the replacement m^ — )■ m(/Jr) (see 
Eq. (lA.ll) ) in the expression ioxo^^yi^^ttx in Eq. ([T]). The option MS_MASS can be applied together 
with options LO, NLO and NNLO. 

Let us briefly discuss the main motivation for this option. So far the mass used in all formulae is 
given as the so-called on-shell or pole-mass which is defined as the location of pole of the quark 
propagator calculated order-by-order in perturbation theory. It is well known that the pole-mass 
is not a well defined concept in QCD since quarks do not appear as asymptotic states in the 
quantum field theoretical description of the strong interaction owing to confinement. In other 
words, the quark propagator does not have a pole in full QCD. A more quantitative analysis 
leads to the conclusion that the pole-mass is intrinsically uncertain of the order of Aqcd (see 
e-g- |[24ll25l ). Since in perturbation theory the pole-mass can be expressed in terms of a short 
distance mass like for example the running mass which is free from non-perturbative effects it is 
advantageous to translate the cross section predictions from the on-shell scheme to the MS mass 
scheme. As a benefit, the convergence of the perturbative series is significantly improved when 
the running mass is used and the extracted numerical value of the top-quark mass is very stable 
under higher order corrections. These observations were first pointed out in Ref. lfT2l . 
In the Hathor program the conversion (5h^h2^ttx{S,mt) — )■ (5hyh2^ttx{S,rn{^r)) has been realized 
now in an easy way allowing a direct evaluation of the cross section using the running mass (see 
also [il8il ). All details are deferred to Appendix lAl 



Option PDF_SCAN 

The option PDF_SCAN allows for the automated computation of PDF uncertainties. In the default 
setting of the Hathor package the PDFs can be accessed with the LHAPDF library Il26ll271l . A 
prerequisite for this option is, of course, that the respective PDF provides a set of error PDFs. 
There are, however, different conventions with respect to PDF uncertainties. 
For instance, there exists the convention of asymmetric uncertainties, a choice adopted by e.g. 
MSTW im and CTEQ jm. Here the error PDFs come in updf pairs (a^,+,a^_), where the 
first element of the pair describes the error of the corresponding parameter in the '+' -direction, 
the second the one in the ' — '-direction. Then, for a given PDF set with a central fit resulting in 
a cross section Oq the systematic uncertainty Ao± is estimated by (see e.g. Il30l ). 



Ao± = - / Y, (max(0,±a;t,+ TOo,±0;t,-TOo))2. (14) 



k=\,iipof 



Eq. (fT4l) is the default of the Hathor package when using the option PDF_SCAN. Following stan- 
dard conventions the PDF uncertainty should be linearly added to the theoretical uncertainty 
from scale variations (parameterizing uncalculated higher orders). 

Other PDF sets, e.g. ABKM f,31il . employ the convention of symmetric uncertainties, where the 
npDF elements each describe the (symmetric) '±'-variation. In this case, the quadratic uncer- 
tainty Ao± is obtained by adding the individual errors quadratically in the standard manner. 



^^± = ./ L (^fc-cfo)^ (15) 

y k=l,npDF 

and the option PDF_SCAN has to be combined with the additional option PDF_SYM_ERR. 
Finally, there exist PDF sets, e.g. B21I33II which simply return a number npof of best fits, where 
typically npop — O(IOO). Then, the PDF uncertainty of the cross section o is estimated by com- 
puting it with each of the best fits and taking the standard statistical average. In such cases, the 
option PDF_SCAN cannot be used for an automated computation of the PDF uncertainty. Within 
Hathor, it of course, always possible for the user to provide own code for the evaluation of the 
PDF uncertainy. 

If, however, a PDF set provides different values of the strong coupling a^ for different error 
PDFs, this is automatically taken into account. 



Before continuing with the description of the Hathor program, let us mention that the package 
offers the possibility for several extensions in the future. With an experimental precision of 5% 
as envisaged at LHC the electro-weak radiative corrections at NLO ifTSUTTl have to be taken into 
account. While for the Tevatron they turn out to be small (less than 1%) they are of the order of 
2% at LHC with a slight dependence on the Higgs mass. Electro- weak NLO corrections can be 
included in a similar manner as the higher order QCD corrections, i.e., with the help of accurate 
fits. 

Also the treatment of QCD radiative corrections leaves room for improvement, e.g. by incorpo- 
rating bound state effects from the resummation of Coulomb type corrections lfT3l[T4| . Finally, 
the design of the Hathor package can in principle also accommodate related approaches for the 
computation of the total top-quark pair cross section beyond NLO, for instance those based on 
soft gluon enhancement in differential kinematics If34l[35]| (see Ref. [[36]| for earlier work). Such 
extrapolations of large logarithms from soft gluon emission in a differential variable (e.g. the 
top-quark pair-invariant mass) to the full partonic phase space introduce systematic uncertainties 
and require a detailed comparison to an inclusive approach such as in Eqs. (fT0l)-(fT2l). 

3 Installation 

In the default setting the Hathor package is based on the LHAPDF library Il26ll to access the 
PDFs. The Hathor package has been tested with the most recent version lhapdf-5.8.3 which can 
be obtained from |http: //projects .hepforge.org/lhapdf| Please follow the instructions 
in the LHAPDF package to install the LHAPDF library. To build and use Hathor first unpack 
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the package using tar xvfz Hathor .tar .gz at the location where one wants to install the 

package. This will create a directory Hathor-1.0. Please create a symbolic link Ihapdf inside 

this directory to the location of ones LHAPDF installation. The contents of Ihapdf should 

contain the LHAPDF installation with the directories: bin include lib lib64 share. Then 

use make to build the Hathor package. Make will build the Hathor (static) library libHathor . a 

which can be used in other applications. In addition an example program main . exe is built. For 

details concerning the example program we refer to Section [51 Note that the compilation is done 

using the GNU compilers gcc, gfortran and g++. The package is known to work also with the 

Intel compiler. In fact we recommend the usage of the Intel compiler since this results in a much 

better performance. However, given that it is not everywhere available Hathor uses the GNU 

compiler by default. [^ 

To run the example one has to tell the system where the dynamic libraries for LHAPDF can be 

found. This is conveniently done using the environment variable LD_LIBRARY_PATH. Using the 

C-shell the statement would be: 

setenv LD_LIBRARY_PATH <path_to_lhapdf_installation>/lib/ 

In addition, one probably needs to specify the location where the grid files for the LHAPDF 

library are stored. Again using the C-shell, the statement would be of the form: 

setenv LHAPATH <path_t o_lhapdf_installation>/ Ihapdf /share/ Ihapdf /PDFsets 

For a detailed description concerning the paths required by the LHAPDF library, we refer to the 

LHAPDF manual. As concluding remarks with respect to the PDF libraries we would like to 

point out that Hathor does not try to handle errors from the LHAPDF library. This is not possible 

since LHAPDF does not provide a detailed error handling. Also note that since LHAPDF is not 

thread safe the same is true for Hathor. 

Hathor is also equipped with a graphical user interface (GUI) written in Java. The GUI makes 
use of dynamic libraries to access the Hathor library within Java, a The dynamic library is 
built with the command Create JavaGui .csh which is included in the Hathor package. The 
command Create JavaGui . csh creates the dynamic library libHathor4 Java.so and writes 
the executable file xhathor which is used to start the GUI. Please note that xhathor sets up 
various paths: i.e. the environment variable LHAPATH is set to . /Ihapdf /share/lhapdf if it 
has not been set yet. In addition the dynamic libraries libLHAPDF . so and libgf ortran . so on 
which the Hathor library relies are preloaded. This step is platform dependent and not alway easy 
to achieve in a universal way. If the graphical interface does not start with xhathor the correct 
setting of the paths is a likely source of potential errors. In that case we recommend to set the 
necessary paths directly in xhathor or to consult the authors for support. 



'For the Intel compiler the user has to adapt the makefile. 

^The same technology can be used to access the Hathor library from Mathematica or Maple. The authors may 
provide the respective interfaces in a future update on demand. 



4 Description 

The entire cross section is calculated inside the class Hathor. This is done in order to avoid any 
possible problem with names used in already existing codes. Inside this class, a two dimensional 
numerical integration is performed in which the PDFs are convoluted with the hard scattering 
cross section. As a numerical integration procedure we use the Vegas algorithm Il37]| .rl Since 
Vegas is a Monte Carlo integration we need to provide random numbers. Those are obtained by 
using the Ranlux algorithm |[38l and we use the implementation available from Ref. Il39l . 
The Hathor class takes as constructor argument a reference to the PDF which should be used in 
the current cross section calculation. In the following we list the publically available function 
calls and option choices together with a short description. 

• Hathor (Pdf & pdf) 

Constructor to build one instance of the Hathor class. The argument is an instance of the 
PDF. In case that LHAPDF is used the corresponding definition would be: 
Lhapdf pdf ("MSTW2008nnlo68cl") ; 
tousetheMSTW2008nnlo68cl set. 

At first sight it might appear strange that we use an additional "wrapper class" as interface 
to LHAPDF. The idea behind this is to give the user the possibility to become independent 
from LHAPDF. This is achieved by inherting the class Lhapdf from the base class Pdf. 
By inhering its own class from the base class the user can thus easily implement its own 
wrapper to whatever PDF set he wants to use. As an example the class MSTW has been 
supplied, which gives direct access to the MSTW set ll28l . (Note that in the MSTW case the 
a^ value is set to 1 since the library does not provide a function to evaluate it. The user has 
thus provide its own a^.) We have observed that in some cases the original code provided 
with the PDF sets is faster than what is provided by LHAPDF. Since the evaluation of the 
PDFs represents a significant part of the calculation the usage of the original PDFs may 
speed up the entire calculation significantly. 

• void printOptions 

Use this function to print the options currently selected via the routine set Scheme ( ) ; 

• void setScheme (unsigned int newscheme) 

Sets the specific scheme in particular perturbative order and renormalization schemes. Pos- 
sible options are: 

- Hathor : : LO to switch on the leading-order contribution. 

- Hathor : : NLO to switch on the individual NLO contribution. 

- Hathor : : NNLO to switch on the individual NNLO contribution (see Section[2l). 

Please also note that for the computation of the total cross section up to e.g. NNLO accu- 
racy, it is necessary to combine the options as in 



^Hathor uses Vegas code which is a C port of the original fortran version 1371 . 
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setScheme(Hathor: :L0 | Hathor::NLO | Hathor : :NNLO) ; 
Other possible options are: 

- Hathor : : MS_MASS to use the mass renormalized in the MS scheme (see Section[2l). 

- Hathor : : LOG_ONLY to keep only the logarithmically enhanced terms (see Section[2l). 
Please note that this option requires also Hathor : : NNLO to be set. 

- Hathor : : PDF_SCAN to switch on the evaluation of the PDF uncertainties. That is to 
say, the error PDFs are also integrated along with the central value. To save comput- 
ing time one may set the accuracy with XS . setPrecision (Hathor: :LOW) to LOW in 
this case. Care has to be taken, though, by the user to ensure that the default PDF 
uncertainty estimate as implemented in Hathor (asymmetric error) complies with the 
conventions of the respective PDF set, as e.g. some PDF sets provide only a symmet- 
ric error. In this case, the additional option Hathor: :PDF_SYM_ERR needs to be set. 
See also the discussion in Section |2l 

- Hathor : : PDF_SYM_ERR invokes symmetric PDF error, if foreseen by the convention 
of the PDF set. 

Please note, that these options can again be combined, e.g. as in 

setScheme (Hathor : :L0 | Hathor: :NLO | Hathor : :MS_MASS) ; 

void setColliderType( COLLIDERTYPE type) 

Sets the hadronic initial state. Use Hathor: :PP to select proton-proton collisions and 
Hathor: :PPBAR for proton-anti-proton collisions. The collider energies are set to the 
default values: 7 TeV in case of proton-proton collisions and 1960 GeV in case of proton- 
anti-proton collisions. If this is inappropriate the values can be changed using the com- 
mand void setSqrtShad (double ecms) , where the center of mass energy is provided 
in GeV. 

void setSqrtShad (double ecms) 
Sets the center-of-mass energy in GeV. 

void setNf(int nf) 

Sets the number of massless flavors to nf. For top-quark physics the default setting is 
yif = 5 and should not be changed. This function may be used when the cross section for 
a hypothetical heavy quark of a fourth family is computed, as Hathor includes the full nf 
dependence of the hard scattering cross section, i.e. it features the formulae for general nf. 
However, please note that the PDFs usually provide a^ in the nf = 5 flavor scheme. So the 
user should be careful with this option (and the interpretation of the results). 

void setCqq (double tmp) 

Can be used to set the constant defined in Eq. (flOl) to a specific value (see Section O. The 

default is Cqq = 0. 



11 



void setCgg (double tmp) 

Can be used to set the constant defined in Eq. (fT2l) to a specific value (see Section O. The 

default is Cgg = 

void setPrecision (int n) 

Can be used to define the accuracy of the numerical integration performed by the Hathor 
package. It sets the number of function evaluations used in the Monte Carlo integration. 
In principle, the user can provide any reasonable integer value. 

Pre-defined values are: Hathor: : LOW, Hathor: : MEDIUM, Hathor : : HIGH. We recommend 
LOW for the PDF scan and MEDIUM for the central value. This should be sufficient for most 
applications. Please note that Hathor : : LOW should give already an accuracy at the percent 
level. For detailed comparisons of theoretical results HIGH may be used. 

double getAlphas (double mur) 

Returns the QCD coupling constant at the renormalization scale mur as provided by the 

(central) PDF. 

double getXsection (double m, double mur, double muf) 

This function starts the cross section calculation for a given top-quark mass and the fac- 

torization/renormalization scales provided as arguments. Unless a specific scheme is set 

through setScheme the default setting is used: 

Hathor: :L0 | Hathor: :NLO | Hathor: :NNLO 

The cross section for the central PDF is returned. More information can be obtained 

through getResult. 

void getResult (int pdfset, double & integral, double & err) 
This function is used to obtain additional information after the cross section has been calcu- 
lated for a specific setting of the renormalization/factorization scale using getXsection. 
The integer value pdfset specifies the respective PDF: for the central value, and 1 to 
getPdfNumber ( ) for the respective error PDF. The result for the central value and the nu- 
merical error from the numerical integration are returned through integral and err. Note 
that err should always be negligible. If this is not the case the precision of the numerical 
integration should be increased through setPrecision. 

void getPdfErr (double & up, double \& down) 

This function returns the PDF uncertainty, if the option Hathor : : PDF_SCAN has been used. 
By default, Hathor assumes an asymmetric PDF error convention. In case of a symmetric 
one, the option PDF_SCAN has to be combined with the option PDF_SYM_ERR (see Section|2] 
and the discussion above). 

int getPdfNumber 

Returns the number of error PDFs currently in use. If is returned the option PDF_SCAN is 

not switched on or the PDF set does not support error PDFs. E.g., in case of the PDF set 

mstw2008nnlo. 68cl the return value would be 40. 
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• void sethc2 (double) 
Can be used to change the value for (hc)^ which is used by Hathor to convert the cross 
sections from GeV ^ to pico barn. The default used by Hathor is: 

0.389379323e+9. 

5 Usage and examples 

A concrete instance of the Hathor class is built using: 

Lhapdf pdf ("MSTW2008nnlo68cl") ; 
Hathor XS(pdf) ; 

The evaluation of the cross section (using the default setting) is then done using 

XS.getXsection(171.,171.*2,171./2) ; 

where the mass has been set to 171 GeV and ^^ = 2 x 171 GeV and /jf = 171/2 GeV. The result 
of the evaluation is obtained through 

XS .get Re suit (0, val, err, chi2a) ; 

Note that XS.getXsection (171.,171.*2,171./2); triggers the numerical integration of the 
cross section. It has to be called first before XS .getResult (0, val, err, chi2a) ; can be used. 
In a typical application we may want to use a lower statistic in the Monte Carlo integration for 
the evaluation of the PDF uncertainty. This is achieved by the following code: 

unsigned int scheme = Hathor: :L0 | Hathor: :NLO | Hathor : :NNLO; 
double mt = 171., muf=171 . ,mur=171 . ; 
double val, err, chi2a, up, down; 

Lhapdf pdf ("MSTW2008nnlo68cl") ; 
Hathor XS (pdf) ; 

XS.setPrecision (Hathor: :MEDIUM) ; 
XS . getXsection (mt, mur, muf ) ; 
XS .getResult (0, val, err, chi2a) ; 

XS.setScheme (scheme | Hathor : :PDF_SCAN) ; 
XS . setPrecision (Hathor : : LOW) ; 
XS .getXsection (mt,mur,muf ) ; 
XS .getPdfErr (up, down) ; 

The central value is calculated with the precision set to MEDIUM. The estimate of the PDF uncer- 
tainty is then obtained with a lower accuracy. 

An example of the usage of Hathor illustrates the calculation with running a mass. It reproduces 
the central curve (NNLO) of the right plot in Figure 5 of Ref. lfT2l . 
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double val, err, chi; 

Lhapdf Ihapdf ("MSTW2008nnlo68cl") ; 
Hathor XS (lhapdf) ; 



XS.setColliderType (Hathor: :PPBAR) ; 
XS.setScheme (Hathor: :L0 | Hathor: :NLO 
XS . setPrecision (Hathor : : LOW) ; 



Hathor: :NNLO I Hathor:: MS MASS 



for (double mt = 140; mt < 181.; mt++ ){ 

XS . getXsection (mt, mt,mt) ; 

XS .getResult (0, val, err, chi) ; 

cout << mt << " " << XS.getAlphas (mt) <<" "<< val << " " << err << endl; 
} 

The typical runtimes for these examples range between seconds and a few minutes and also de- 
pend on the chosen compilers. E.g. on an Intel 3.00GHz QuadCore PC we have obtained with the 
options NNLO, PDF_SCAN (PDF set MSTW2008nnlo68cl) and XS . setPrecision (Hathor : : LOW) 
the result for the cross section after 64 seconds using the gfortran compiler, and, 41 seconds 
respectively, using Intel's ifort compiler. 

The Java GUI is invoked by the command xhathor (see the discussion in Sec.[3]for the installa- 
tion). A screenshot is displayed in Fig.[T]and the input is self-explanatory with the help of Sec.|4] 
for the description of all options. 
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Figure 1 : Screen shot of the Java graphical user interface for Hathor. 
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6 Conclusions 

Top-quark production at hadron colliders is at the edge of becoming a precision measurement 
requiring accurate theory predictions. Hathor is a fast and flexible program for the computation 
of the total cross section of hadronic heavy-quark pair-production. It takes into account the 
latest theoretical developments through a variety of options, it allows for separate variations of 
all scales and can be used with a large number of PDFs through the LHAPDF interface. As a 
novelty, Hathor offers predictions in different renormalization schemes for the heavy-quark mass 
(pole and MS) and it can also be applied to a hypothetical fourth quark family assuming standard 
QCD couplings. Thus, with these functionalities Hathor can serve as a reference for future cross 
section calculations. 

Hathor typically runs in a few seconds up to a few minutes (depending on the chosen options, 
e.g. extensive PDF scans) on standard desktop or notebook PCs. The Hathor package can either 
be used as a stand alone program or, as a small library, it can be easily integrated into existing 
code, e.g. for experimental analyses. 

Hathor is publicly available for download from [|40ll or from the authors upon request. 
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A Total cross section with a running mass 

The starting point is the relation between the on-shell mass and the MS mass: 

m, = m(A/r)(l+fl4+«?2) , (A.l) 

with Qs = ds/n. If the decoupling a^"/=^ — )■ aj'^f^^ is performed at m(/7r) the coefficients are 
given by 

1 = ^+^, (A.2) 

307 2 1 509 47 ^ 

2 = |^+2C(2) + iC(2)ln2-ic(3) + 9^« + ^<^ (A.3) 



(^ + 5«2) + i^ + n'^)"/ + 5EA(".,M) 



which are known from Refs. [|4TI - |43ll and I = ln(/j^^/m(^^) ). Note that the coefficients di de- 
pend in general on the renormalization scale. Using m{m) instead of rn{iJ.r) the formulae sim- 
plify significantly. The full renormalization scale dependence can be restored at the end using 
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renormalization group arguments. A{mi/mt) accounts for all massive quarks m, lighter than the 
top-quark. For all light quarks we set mi = so the sum in Eq. (IA.3I) vanishes. 
To convert the cross section to the MS mass scheme we start from the hadronic cross section 
expanded in a^: 

Expressing m; in terms of m{Tn) and expanding in a^ we obtain 



(A.4) 
(A.5) 
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The derivatives of the LO cross sections with respect to the mass can be written in the following 
form: 



Jo(0)( 
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i{m) J Amimf"" y ds '■'"''' ""-"y dmt 

mi =m(jn) 

where a summation over the contributing parton channels is understood. Note that in Eqs. (IA.5I) . (IA.6I) 
and (IA.7I) the renormalization scale is set to ^r = m{m). The required derivatives of the LO par- 
tonic cross sections with respect to the mass are easily obtained from Eqs. ^, ([5]): 
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For the first derivative of a'^^ we obtain a similar result: 
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(A.8) 
• (A.9) 

(A.10) 
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with 



(1) / „ ^,^. ^f 



~(1). — /— NN 3 (1)/ 

a}j'{s,m{m)) =ijf—a]j' I 5,m(m) 



m(m] 



Using Eq. ^ the contribution d| ■ (5, m{m)) can be written as 



2~(1) o «v ("^H) ^(ll) 

^ c O; ; — 8 ^ J- J 



IJ 



m{m) 



Since the luminosities are only known numerically the derivatives are evaluated using 

—L,j{s,S,Hf) = — {L,j{s + d,S,^f)-L,j{s-d,S,^f))+0{d^). 



(A.ll) 



(A. 12) 



(A.13) 



The results presented so far are only valid for /Ur = m{m) . Using 

as{m{m)) = a,{pr) (l +4Ti\{^ir)Ljiq + {4%Yasipr)HlLR + Q^Ll)) , (A.14) 

with L^ = \n{fjr^ /Tn{Tn) ) it is easy to restore the complete renormalization scale dependence. 

B Scaling functions 



Here we present the expressions for the scaling functions as implemented in the program Hathor. 
At Bom level, 
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where p = -^/l — p and p = 4-mf/s. At NLO the functions /^ have been described through 
precise parametrizations with per mil accuracy and the following ansatz ||T2]| : 
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where nf denotes the number of light quarks and the complete n /^-dependence has been kept 
manifest. The constants c/^ read af = 0.03294734 and aff = 0.01875287 and the fit functions 
/z([3,ai,...,ai7) andhgg{^,ai,...,ai5) are given in Eqs. (IB. 181) . ( IB. 191) together with a list of all 
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parameters in Tabs.[I]-[3l The exact expressions for scale dependent functions f^j have already 
compact analytical form containing at most dilogarithms. They read |[T1 
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At NNLO the functions f>. ' and /j- ■ are known exactly lfT2]| . The fits to the scaling functions 
generally have per mil accuracy with exceptions in regions close to zero, where an accuracy 
better than one percent is kept. 
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Table 1 : Coefficients for fits of the qq scaling functions. 
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Table 2: Coefficients for fits of the gq scaling functions. 
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.(10) 


(21 

Jss 


) 


A22 
Jgg 


) 


i 


at 


bi 


Ci 


bi 


Ci 


1 


-8.92563222 


-4.18931464 


0.12306772 


0.01222783 - 


-0.00380386 


2 


149.90572830 


82.35066406 


-2.75808806 


-0.77856184 


0.08757766 


3 


-140.55601420 


-87.87311969 


3.19739272 


1.33955698 - 


-0.10742267 


4 


-0.34115615 


9.80259328 


-0.56233045 


-0.59108409 


0.02382706 


5 


-2.41049833 


-1.12268550 


0.03240048 


0.00248333 - 


-0.00099760 


6 


54.73381889 


29.51830225 


-0.92541788 


-0.23827213 


0.02932941 


7 


90.91548015 


48.36110694 


-1.57154712 


-0.38868910 


0.04906147 


8 


-4.88401008 


-7.06261770 


0.35109760 


0.28342153 - 


-0.01373734 


9 


-0.17466779 


-0.08025226 


0.00227936 


0.00010876 - 


-0.00006986 


10 


13.47033628 


7.01493779 


-0.21030153 


-0.03383862 


0.00658371 


11 


22.66482710 


15.00588140 


-0.63688407 


-0.29071016 


0.02089321 


12 


4.60726682 


3.84142441 


-0.12959776 


-0.11473654 


0.00495414 


13 


-67.62342328 


-47.02161789 


1.91690216 


0.98929369 - 


-0.06553459 


14 


-9.70391427 


-8.05583379 


0.26755747 


0.24899069 - 


-0.01046635 


15 


65.08050888 


47.02740535 


-1.86154423 


-1.06096321 


0.06559130 


16 


5.09663260 


4.21438052 


-0.13795865 


-0.13425338 


0.00551218 


17 


-20.12225341 


-14.99599732 


0.58155056 


0.35935660 - 


-0.02095059 



Table 3 : Coefficients for fits of the gg scaling functions. 
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